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Abstract 

We study an ultra-cold atomic Fermi Gas with the atom-atom interaction modulated periodically 
in space. A novel ground state with cooper pairs occupying non-zero center of mass momenta 
is found. Such a state is closely related to the state proposed by Fulde, Ferrell, Larkin, and 
Ovchinnikov(FFLO). The resultant single particle excitations with momenta along the direction of 
the modulation shows multiple-gap structures. Such a system can be realized in experiment with 
a spatially modulated Feshbach resonance. Experimental signatures of such a state are discussed. 
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In several recent experiments, polarized Fermi gases have been realized^]. These pro- 
gresses drive a new wave of studies of the so called FFLO state ^0]. Such a state is of great 
interest because of its implications in many systems, including color superconducting quark 
matter, dense nuclear matter, and superconducting systems in solid states^]. The observa- 
tion of such a state, however, is still elusive because it exists only in a very narrow parameter 
regime in cold atomic gases. In addition, it is not completely clear whether such a state 
favors occupation of multiple or singlej^] center of mass (COM) states. In this paper, we 
propose to study a ground state closely related to the FFLO state. At zero temperature in 
a homogeneous system, it is the only ground state of the system and the study of such a 
state should shed light on the study of the FFLO state. 

We consider a fermion gas with the atom-atom interaction varying in space. For an 
ultra cold atomic gas, the effective atom-atom interaction is, V (p) = ATragfi^S (p) /m where 
p is relative coordinate between two particles, m is the mass of the atom, and is the 
s-wave scattering length. Near a Feshbach resonance, the scattering length can be described 
by 0] aeff = ttbg — ^^2y^ Ifl'ol^j where a^g is the background scattering length, and the 
second term is the contribution from the Feshbach resonance nearby. The constants Qq and 
Vq are the coupling strength and energy detuning between the scatte ring channel and the 
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can be controlled by an external magnetic field. In an optical Feshbach resonance 

HQ 

can be tuned by laser frequencies, and the free bound coupling can be tuned by the laser 
intensity. A spatially dependent scattering length can be achieved by applying spatially 
varying external fields. This gives an effective interaction, y(r,p) = 47ras (r) (p) /m, 
that depends on both the relative coordinate and the center of mass (COM) coordinate r. 
This direct substitution of the scattering length locally requires that the scattering length to 
vary only in a scale much larger than the local scattering length. In particular, we consider 
the simplest form of such a interaction, 

V(r,p) = [(7o + (7cos(ko-r)]5(p) (1) 

The constant qq term, with qq < g < is added to guarantee that the system stays in the 
BCS side of the resonance. Such a system, which to the best of our knowledge have not 
been studied before, is only feasible with recent advances in atomic physics. A relatively 
obvious phenomenon that arises from such a spatially modulated interaction would be the 
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modulation of the atomic density. Atoms are prone to congregate in locations with maximum 
attractions. A much more interesting phenomenon resides in the pairing of atoms. In a 
typical BCS theory, the ground state wave function consists of Cooper pairs with their 
COM momenta equal to zero. The interaction we propose (Q) creates coherences between 
pairs with COM momenta differing by ko in order to lower the free energy of the system. The 
ground state thus includes Cooper pairs with zero, as well as ±ko, ±2ko COM momenta. 
However, we find that the probability of occupying higher COM momenta is low due to 
disadvantage in kinetic energies and the constant attractive interaction qq. Our ground 
state is closely related to the FFLO state discussed in the literature [3] In the FFLO 
state, due to the mismatch of the chemical potential for spin up and spin down particles, 
particles can pair with non zero COM momenta to lower the energy. Before proceeding with 
our analysis, we emphasize that although our ground state is similar to the FFLO one, the 
spatial symmetry is spontaneously broken in the FFLO case, while in our case the spatial 
symmetry is broken by an external field. 

We work in a grand canonical ensemble and use a single channel Hamiltonian which is 
appropriate for the case of a broad Feshbach resonances. 



where /z is the chemical potential, and ipa (r), ^1 (r) are the fermionic annihilation and 
creation operators with spin index a at position r. In the interaction part, the relative 
coordinate is already integrated out. The V (r) is the standard renormalized interaction^] 
from our bare interaction (^. We solve this system in a mean field level at zero temperature 
using the Bogoliubov de Gennes(BdG) transformation [l^ 





{sk - /u) Mk + Uk'Uk-k' + Ak'Wk-k' = Evk 
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A summation convention is used for any momentum that appears twice on the left. 
We take the gap as A/j^ (r) = — (7 (r) (■?/'^ (r) ^/'q (r)) and Hatree potential as f/aa (r) = 
g (r) i^l)]^ (i") V^Q (r)). The spin index can be neglected with no confusion due to the sym- 
metry between spin up and down. In the following we take kg along z direction. With this 
choice, A (r) , U (r) are functions of z only. In momentum space, we can solve the single 
particle problem in a set of subspaces corresponding to constant kx and ky. The Fourier 
transforms of u{y) ,v (r) , (r) A (r) and U (r) are defined as Mk, fk, S'k.Ak and f/k respec- 
tively. We can take m (r) , f (r) and Akto be real without loss of generality. Since we are 
dealing with a problem invariant under reflection r — r, we can also take Ok, 6k, Ak and 
f/k to be real, ^k = h'^k'^/2m is the kinetic energy of a particle with wave vector k, and V is 
the quantization volume. We note that near a Feshbach resonance the chemical potential, fi, 
need to be self-consistently determined by an average density of the Fermi gas n in equation 

We solve the BdG equations self-consistently by numerical iteration. The self-consistent 
gap is found to have only three components in momentum space, namely, kz = 0, ±ko. 
Due to the form of interaction we chose (|T|). the zero momentum component gap is always 
larger than the component of momenta ifco and therefore the gap is nonzero everywhere 
in real space. This is different from the pairing gap discussed in the FFLO papers 
where the gap has only k^ or i/cg components. The pairing of the many body ground state 
is given by (ctkT'^q-ki)- Note that in the k^ — ky plane, the pairing occurs with opposite 
momenta. In the z direction, however, there are non zero COM q values. We study the 
pairing of atoms with their k^ = ky = without loss of generality. It can be shown that 
(^fcztS-fczi) ~ the numerical result is presented in figure(IT^). 

We find that only states with COM momenta q = 0, ±ko are occupied. The higher COM 
momentum states with integer multiples of kg have negligible wight. Given the magnitude 
of the coupling between different COM states components, it is surprising that atom pairs 
do not occupy states with COM momenta being higher harmonics of k^. Actually this result 
can be explained by the nonlinear property of the BdG equations. Atom pairs tend to 
occupy as few COM states as possible because the interaction energy is proportional to the 
square of {(ik^-\Cig-k^i) ■ By occupying only these three COM states, the pairing energy is 
maximized in magnitude. The pairing occurs mostly for atoms near the Fermi surface. For 
a particular kz, the pairing amplitudes for the COM q = ±kQ are generally not equivalent. 
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FIG. 1: Pairing amplitudes ak^iCbq-kzl kx = ky = 0. Plot (a) corresponds to the ground state, 
(b) to the first excited state, (c) to states right before the second gap, and (d) tp the state right 
after the second gap. Only momenta close to the fermi surface are shown. The arrows are used as 
guide to the eye to pair amplitudes where significant changes taken place. We have used a units 
with fcrmi momentum and fcrmi energy to be one. The numerical values for the parameters are: 
go = 20, g = 15, and ko = 0.2. The self-consistant gaps and chemical potential are found to be, 
Ao = .13, A±ko = -06 and /x = 0.85. 

i.e. (afc^Tflfco-fca) 7^ (^fe^T'^-fco-fca) • This is because for ^ 0, either k^ - k-^ or -k^ - k^ 
is closer to the Fermi surface. Since pairs can occupy several COM states coherently, one 
particular atom, say aj^p can form pairs simultaneously with several particles, say a_kj., 
flko-kj.,and a_ko-kj,- Right at the Fermi surface, pairings with opposite momenta dominate. 
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and therefore the amphtudes (akj^T'^ko-k^^j,) and (ok^t'^-ko-kf i) are small. Slightly away 
from the Fermi surface, the pairing with opposite COM momenta is not as strong and one 
starts to see pairings with nonzero COM. This feature of the ground state pairing is closely 
related to the single particle excitation spectrum discussed below. 

The eigenenergies that we found by solving Q (jH) are the single particle excitation spec- 



trum of the many 



3ody state. The excitation is always gapped as opposed to the state 



studied by FFLOj^j^. The quantity A (r) is positive everywhere in real space with its 
minimum at Aq — A^,,. This minimum is the lower bound of the excitation gap. Above the 
gap, the excitation spectrum is continuous due to the continuous distribution of kx and ky. 
For a particular set of k^ and ky, for example k^ = ky = 0, additional gaps emerge due to 
the periodic modulation of the interaction along z direction [see Fig(j2))]. 

The modulation therefore induces band gap structures in quasi-particle excitations. More 
interestingly, we can understand the excitation spectrum in a pair breaking picture. The 
pairing amplitudes in an excited state \n) are, {n |afc|aq_fc|| n) = (0 \ak^aq-ki \ 0) — u^k '''^k^-q- 
The pairing amplitudes for relevant states are shown in Fig Ei)- The pair breaking 

happens simultaneously around both —kp and kp- The energy lowered by this superposition 
is a small effect, and in the following we focus on one side of the Fermi surface. For the 
lowest excited state the change in the pairing amplitudes occurs right at the Fermi 
surface. The amplitude ak^jCt-kFi is almost completely destroyed in state One the 
other hand, recall that in the ground state, there are no pairing amplitudes for ctkFT'^ko-kFi 
and akpta_ko-kp.j,, due to the dominant pairing of a^pi^a^kpi. Now that the ak^ta_k^| pair is 
broken, the particle a|.^^ is free to pair with particles a±ko-kj^i ^^is decrease the energy 
by a amount of A^g Since the excitation is close to the Fermi surface, the kinetic energy 
contribution can be neglected. Furthermore, the Hartree term is almost the same for the 
ground and excited states and it can be neglected too. Therefore the emergence of gaps can 
be attribute to the pairing energy only, namely 2Ao'U^"''t>["''' + 2AkQU±ko+kVk ■ The gap 
energy can then be estimated to be 0.09, close to the observed value 0.10. The excitation 
energy increases continuously as the excitation gets further and further away from the fermi 
surface. The states with energies right below and above the second gap have sudden flips of 
the signs for the pairing amplitude ak'^a±ko-k'i- Here the k' are the particular wavevectors 
where pairing amplitudes ak'ia±ko-k'i undergo significant changes with respect to the ground 
state. We can repeat the analysis above and estimate the energy cost for the sign flipping 
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FIG. 2: Plot of the single particle excitation energy. Here we only look at the eigen-excitations 
with kx = ky = 0. The n labels the eigenstates. We have used same numerical values for relevent 
parameters as figure (0). 

to be 0.09, close to the observed value 0.11. Therefore, the phase flips of the pairs with 
nonzero COM momenta give rise to the second gap. Similar observations can be made to 
locate gaps at higher energies. However, as the kinetic energy increases, the gap structure 
gets more and more obscure. Actually, with parameters used in the calculation, only two 
gaps are observed. We would like to emphasis that, after summation over the contribution 
from different and ky, the additional gap we discussed above, is not a true second gap 
in the excitation spectrum. However, it nevertheless induces a sharp change of density of 
states for energies close to the gap. Such effects, arising from the rich structures of the pair 
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wavefunction, can be observed in experiments. 

The parameters we use in the calculation are based on recent ^Li experiments QlPl- We 
take the fermi energy Ep = {Stt^u)^^^ ~ 3^k and momentum kp to be the unit energy 
and wave vector, respectively. In this units, the background scattering length a^^ ~ —0.5 
and the modulation of the interaction (Q) /cq ~ -2. Here we estimate to be the order of 
laser wavevector. It is smaller than kp and this choice automatically satisfies the condition 
ko ^ |27r/a|, necessary for the validity of equation (P). The constant coupling strength is 
about qq ~ —20. The coupling between different COM g is chosen to be 15, smaller in 
magnitude than qq. This coupling can be realized in experiment both by external magnetic 
and optical fields. However, it is not easy to generate spatial variations of the magnetic 
field in a scale as small as the sample size. Optical Feshbach resonance, as demonstrated 
in the experiment is a more promising way to generate the spatial modulation of the 
interaction. By coupling the incident channel of two atom scattering to the molecular excited 
state close channel ( 111] . one can induce the periodical spatial modulation by an optical lattice. 
The laser field frequency needs to be chosen to be far away from the resonance of a single 
atom to minimize the induced single particle effects. 

The observation of this novel pairing in the ground state involves measurements of the 
COM momenta of pairs. This can be done by pair-wise projecting Cooper pairs into 
molecules by a fast sweep from BCS side of the resonance to the BEC side. After this pro- 
cedure, the molecules of the condensate should coherently occupy momentum states q = 0, 
q = ±ko, which can be measured by a time of flight image of molecules. Such measurements 
have already been used in several BEC-BCS cross over experimentsl^l. Another way of ob- 
serving the pairing correlation, recently demonstrated experimentally |l5i] . is to measure the 
shot noise correlations. For the excited states properties, one needs to measure the radio 

and recently used to measure 

the fermi gas pairing gapj^. In experiments, two of the hyperfine levels of atoms, |1) , |2) 
are identified as the spin up and spin down states. A probing RF field couples one of the 
hyperfine levels, say |2) , to a third hyperfine level |3). The transition rate can be evaluated 
directly from our BdG solution of the quasiparticle excitations. However, Such calculation 
turns out to be numerically challenging. This requires to discretize the momentum space 
in a very small grid. Here, we adopt a Local Density Approximation(LDA) approach. The 
calculation should be valid for ko -C kp, which is exactly the situation we are interested 
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spectrum is calculated locally and the final result is a summation of 
l^, R{uj) = J d^rD {Sk) [Sk (w)] . Note that the local gap and 



in. In LDA, the RI 
local contributions 

chemical potential are given by our self-consistent BdG solutions. Here uj is the detuning 
of the laser frequency from the frequency difference in |2) and |3) , and Q is the coupling 
between them. D [sk) is the free particle density of states and it is proportional to ^/ek. In 
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side the step function, the Sk (cu) = "2^,'''''' + (^) should be greater than zero, which gives 
the threshold of the excitation. The chemical potential of state |3) is taken to be zero and 
the RF spectrum is shown in figure 0- 

We note two features of the spectrum. First, the threshold is shifted to a lower energy 
comparing with the uniform case. This corresponds to a local gap A (r) = Aq — A^^, which 
is consistent with the observation that the lowest excitation states correspond to break a 
COM zero pair and increase the non-zero COM pairing. Second, there is an extra peak at 
the higher energy side of the spectrum, which corresponds to a local gap A (r) = Aq + A^^. 
This is consistent with our analysis that the second gap originates from sigh flips of the 
non-zero COM pairs and the magnitude is proportional to the A^q only. 

In conclusion, we studied a fermi gas with spatially modulated interaction in a mean 
fleld level at zero temperature. We also discussed its experimental realization and detection. 
Such a state has a periodical modulation of the order parameter similar to the FFLO states. 
Even though we considered the spatially varying interaction along only one direction, our 
analysis can be easily generalized to a system where the interaction is modulated in three 
directions. In this case the second gaps along three directions can overlap and produce true 
additional gaps of the quasiparticle excitation. This should produce more pronounce signals 
in the RF spectrum. 

HF would like to thank P. R. Berman for discussions on the modulation of scattering 
length and thank L.-M. Duan and W. Yi for several useful comments. This work is done 
under flnancial support of NSF0244841 and FOCUS Center grant. A.G.R acknowledges 
support from the Research Corporation, Cottrell College Science Award. 
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FIG. 3: (color online) The RF spectrum is plotted for a homogenously interacting system and our 
spatially modulated interacting system. Notice that there is an extra peak in the RF singal for 
the system with a spatially modulated interaction. We have again used the numerical values as 
Figure!^ for various parameters. 
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